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Q^l We explicitly determine the residues and orders of all the ultra-violet 

(UV) and infra-red (IR) renormalon poles in the Borel plane for the QCD 
vacuum polarization function (Adler D-function), to leading order in an ex- 
pansion in the number of quark flavours, Nf. The singularity structure is 
rS \ precisely as anticipated on general grounds. In particular, the leading IR 

c^ \ renormalon is absent, in agreement with operator product expansion ideas. 

There is a curious and unexplained symmetry between the third and higher 
UV and IR renormalon residues. We are able to sum up separately UV 
and IR contributions to obtain closed form results involving (^-functions. We 
argue that the leading UV renormalon should have a more complicated struc- 
ture than conventionally assumed. The disappearance of IR renormalons in 
flavour-saturated SU(A^) QCD is shown to occur for A^ = 3, 6 or 9. 



1 Introduction 

The topic of the large-order behaviour of the expansion coefficients in per- 
turbative field theory calculations continues to be vigorously pursued in the 
literature [1-9]. Notwithstanding the importance of establishing to what ex- 
tent physical observables can be formally reconstructed from their divergent 
perturbative expansions, there is the more pragmatic motivation of the need 
to assess the reliability of the growing number of higher-order perturbative 
calculations in QCD and QED. 

In this paper we wish to make use of recent progress |^, §] in exact all- 
orders QED calculations to leading order in the 1/A^/ expansion, with Nf 
the number of fermions, to explicitly determine the singularity structure 
of the perturbative QCD vacuum polarization function in the Borel plane. 
On very general grounds one anticipates branch point singularities evenly 
spaced along the positive and negative real axis in the Borel variable [§]. 
Those on the positive axis, referred to as infra-red (IR) renormalons, are 
supposedly correlated with the absence from the formal perturbation series 
of infra-red non-perturbative effects, vacuum condensates, present in the op- 
erator product expansion (OPE). They are responsible for fixed-sign factorial 
growth of the series coefficients and represent a genuine ambiguity in recon- 
structing the physical observable from the formal perturbation series. Those 
on the negative axis, so-called ultra-violet (UV) renormalons, correspond to 
alternating-sign factorial growth of the series coefficients and do not prevent 
the reconstruction of the observable by Borel summation. 

Whilst the above singularity structure is well-motivated theoretically, 
there have been various problematic issues. In particular the connection 
with the OPE suggests that the leading IR renormalon should be absent for 
the case of the QCD vacuum polarization function since there is no relevant 
operator of dimension two; the first contribution being the gluon condensate 
of dimension four. This conclusion has been questioned on various grounds 
by several authors [|T^, 0. 

The leading asymptotic growth of the perturbative coefficients will be 
determined by the Borel plane singularity nearest the origin; for the case of 
the QCD vacuum polarization function this is the first UV renormalon. We 
point out that the conventionally expected structure of this singularity, with 
a single branch point exponent, would enable one to obtain the asymptotic 
growth of the coefficients to all orders in the 1/iVj expansion given an exact 



large- Nf result. We suggest that this is unhkely and indicate a more coniph- 
cated structure for the first UV renormalon in accordance with recent results 
of Vainshtein and Zakharov obtained using their "UV renormalon calculus" 

i- 

We shall show that the actual singularity structure of the QCD vacuum 

polarization function is precisely as expected; in particular the leading IR 

renormalon singularity is indeed absent. This has also been noted for the 

singularities in the QED vacuum polarization function in reference p. We 

further demonstrate that there is an unexpected symmetry between the third 

and higher UV and IR singularities. We are able to sum up the UV and IR 

contributions separately to obtain a closed form result involving (^-functions. 

We finally show that in SU(3) QCD, with A^/=15 or 16, the IR renor- 
malon singularities are absent ||T2| , p!3|| ; and that they first vanish when the 
instanton/anti-instanton singularity becomes leading. The requirement that 
this happens for an SU(A^) theory uniquely selects A^ = 3. 

The organisation of the paper is as follows. In section 2 we shall introduce 
the Adler D-function and define its 1/A^/ expansion. We shall also discuss 
the exact large- A^/ QCD result for its coefficients 0. In section 3 we define 
the Borel plane and review the anticipated singularity structure. We discuss 
reasons for which one may expect a more complicated singularity structure 
for the first UV renormalon and make contact with the results of reference [^ . 
In section 4 we shall explicitly determine the Borel plane singularity structure 
for the exact large-iVj result for D; and show that the anticipated UV and 
IR renormalon poles are present. In section 5 we discuss the vanishing of the 
IR renormalons for certain values of Nf in SU(A^) QCD. Section 6 contains 
our conclusions. 



2 The Adler D-function and the l/Nf Expan- 
sion 

We shall be interested in the SU(A^) QCD vacuum polarization function with 
Nf flavours of massless quarks. 



To avoid an unspecified constant we shall actually focus on the related 
Adler D-f unction, 

Dm = -iQ'd^n(g^) (2) 

where Q^ = —q^ is the spacelike Euclidean squared momentum transfer. This 
quantity is related to the experimentally-relevant R-ratio in e"'"e~ annihila- 
tion, 

cr(e"'"e~ -^ hadrons) , . 

-^ = — n^^ + -\ (3) 

(T(e+e -^ n+n ) 

where, taking s to be the physical timelike Minkowski squared momentum 
transfer, R and D are related by the dispersion relation 

In QCD perturbation theory we have 

D{Q') = d{R) E Q? (l + IcfD^ + (e Qf] D (5) 

where Qf denotes the electric charge of the quarks and the summation is 
over the flavours accessible at a given energy. d{R) is the dimension of the 
quark representation of the colour group (here d{R) = N). We define the 
SU(A^) Casimirs Ca = N,Cf = (N^ - l)/2N. 

The correction to the parton model result has the perturbative expansion 

D = a + dio^ + d2a^ + . . . + 4a^+^ + . . . , (6) 

with 'a' the renormalization group (RG) improved coupling Q;s(yU^)/7r. The 

D contribution first enters at O(a^) due to the existence of diagrams of the 
"light-by-light" type. Our interest here is in the asymptotic growth of the dk 
coefficients in large orders. 

The RG-improved coupling alfi"^) will evolve with renormalization scale 
/i^ according to the beta-function equation 

— -^ = -6a2(l + ca + C2a2 + ...) (7) 

dm/x 



where b and c are universal with [0] 
^ _ (llC^-2iV^) 



6 

7Cl 



UCaCp 5 3 ■ 

y-^ + 4^^ + 4^^ 



(8) 



C2 and higher coefficients are renormahzation scheme (RS) dependent. We 
shall usually consider the MS scheme with /i^ = Q^. For the R-ratio there is 
an analogous expansion for the quantity R with perturbative coefficients r^, 
defined as in equations (5) and (6). The dispersion relation (4) means that the 
Tfc are directly related to the dk- For instance ri = di and r2 = d2 — 7r^6^/12. 
The vr^ terms arise due to analytic continuation. Given knowledge of the 
asymptotic growth of the dk one can obtain that of the r^ using equation (4). 
We shall continue to focus on the d^ for the moment. 

The coefficients dk will be polynomials of degree k in Nf-. 



4 = S^^N'f + dt'^Nt' + ... + d 
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/ 



(9) 



where each term is a sum of multinomials in Ca, Cp and Nf of degree k so 
that d^l^ has the structure C^'^~'^Cp (note the prefactor of Cp in equation 
(5)). The first two coefficients di and d2 have been computed |^ and the 
result using the MS scheme with /i^ = Q^, expanded in Nf as in equation 
(9), is 



di 
d2 



11 2 

12 + 3^' 
19 

27 
29 



/41 
Nf + CA[- 



11 



-0 



-'.Cp 



Ca) iVj + Ca 



/151 

V162 
^ , -^ 19 10 

+^^ ' "96 + Y^' ~ y^' 



970 224 5 

"sT ^ yr^^ ^ 9^^ 
/ 90445 



+CaCp 



127 143 



^ ' ^ ^ V 2592 

55 



Nf 
2737 



_55 

lOS" 18 



48 



12 



Cs + ^Cs +^F -777^ 



23\ 



32/ 



(10) 



Here (p denotes the Riemann zeta function, 



1=1 



(11) 



1/A^/ expansions as in equation (9) have been widely used in the past 
in the investigation of large-order behaviour and renormalons 0. As we 
shall emphasise, it is actually more useful for these purposes to consider an 
expansion in powers of b. We can write 

4 = di'¥ + df-%'~' + . . . + 4°^ (12) 

The leading coefficient in the 1/A^/ expansion is exactly related to that in 
the "1/6 expansion" with 

df = i-im'dP . (13) 

As before, we can write out the known di and d2 coefficients now expanded 
according to equation (12): 

V 4 ^V 12 8 

/151 19 \, 2 ^ /31 5^ 5^ 

+c.(|-Hc3 + ioc.). + c5(-I|-C3) 

,C.C.(-?| + HC3).CJ(-|). (14, 

We note in passing that the "1/6 expansion" has a somewhat more compact 
structure than the 1/A^/- In particular the (3 terms in di, which are present 
in all orders of the 1/A^/ expansion, are now present only in the leading 
term in the "1/6 expansion"; and the C5 terms present in all but the leading 
coefficient in the 1/A^/ expansion are now present only in the ^2 coefficient 
in the "1/6 expansion" . 

Continuing progress in applying the 1/iV/ expansion in QED [|l^ has led 
Broadhurst to an elegant generating function for the leading order (large- 
Nf) coefficients of the QED Gell-Mann-Low function (MOM scheme beta- 
function) m]. 

32-" / d 



2 \dx 




x=l 



(15) 



where 

^^ 3(l + a;),t^2('^'-^')' ■ 



\1/|^1 can be explicitly evaluated in closed form ||^ 

n + A 



^i^J [n - 1) 



2n + 4 



{n-2)\ (-3)"-! 

16 



f >s>0 



(17) 



Using this result one can then obtain the leading-order large- A^/ result for 
the QCD Adler D-function. In the MS scheme with /i^ = Q^ one has [^ 

where Tf is a group theory factor; Tf = 1/2 for the standard fermion repre- 
sentation. The (—5/9)"* factors enter since one is converting from the MOM 
scheme Adler function to that in the MS scheme. The results of (18) and (17) 
are in agreement with the exactly known coefficients di and ^2 in equation 
(10). 

Our aim is to make use of the exact large-A^j result of equation (18) to 
obtain as much information as possible about the singularity structure of 
the QCD D-function in the Borel plane; and hence about the large-order 
behaviour of its perturbative coefficients. To this end we shall begin by 
reviewing what can be inferred on very general grounds about this structure; 
and then we shall compare the exact result with these expectations. 



3 The Borel Plane and Renormalons 

Consider a general quantity D, calculated in perturbative field theory with 
coupling a, 

D = {a + dia^ + t/sa^ + . . . + 4^''+^ + . . .) (19) 

Using the integral representation of k\ 

k\= / rfte-y (20) 

Jo 

we can formally write 

'^ rl roo 

D=ya^+^^ dte-'r (21) 



with do = 1. Then exchanging the order of summation and integration 



•^•J m=0 '"• 

= / d.e-^/'^E 



m=0 



m=0 "^• 



/•oo 

= / d2e-"/"5[D](2) (22) 

Jo 

where we have set z = ta. B[D]{z) is called the Borel transform of D{a) and 
is defined by 

B[D]{z) ^Y.-^- (23) 

The idea will be that even if the coefficients dk grow like A;!, as is expected in 
field theory, the Borel transform defined in equation (23) may still be defined 
and hence the integral representation of equation (22), the Borel sum, may 
exist. 

As an example suppose that dk = {—l)'^k\; then B[D]{z) = 1 — z + 
z"^ — z^ + . . . = 1/(1+2;), where we have assumed analytic continuation 
along the whole real line. More generally, if dk = {l/zi)^k'^k\ (7 > 0), then 
B[D]{z) has a singularity proportional to {z — Zi)~'^'~^; so if 7 is a positive 
integer we have a pole; and for non-integer 7 a branch point in the z-plane 
(Borel plane) aX z = Zi. If all the singularities are located off the positive 
2;-axis it may be possible, if certain analyticity properties of D are satisfied, 
to reconstruct D{a) from its formal divergent series using the Borel sum of 
equation (22). This, however, is not our interest here. We want to use the 
z-plane singularities to encode the large-order behaviour of the perturbative 
coefficients. 

Returning to the specific example of the Adler D- function in QCD, let 
us now ask on rather general grounds where the singularities Zi could be 
located. In the large- A^j limit we know that, if the dk have factorial growth, 
asymptotically we must have dk ~ Nfk\{ ). Similarly we can consider a 
large- A^ limit where necessarily dk ~ N k\{ ). This follows since the k^^ order 
Feynman diagrams for dk necessarily have factors which are multinomials in 
Nf, Ca and Cp of degree k. We therefore have singularities at positions 
Zi ~ 1/A^/ in the large- A^/ limit and Zi ~ 1/A^ in the large- A^ limit. If 
singularities are present and visible in both limits then the simplest possibility 



is Zi ~ l/{AN+BNf) involving some unspecified linear combination of A^ and 
Nf. These are the renormalon singularities and in fact they lie at Zk = 2k /b 
where k = ±1, ±2, ±3, .... There are also singularities due to instanton/anti- 
instanton solutions of the classical equations of motion p, |1^ |18| . These lie 
at Zk = 4:k, k = 1, 2, 3, . . .. Since their positions are independent of N and 
Nf, the above arguments suggest that they are invisible in the large- A^ and 
large- A^j limits. In fact they are invisible at all orders of the 1/A^ and 1/iV/ 
expansions, so we shall not learn about them from the exact large- iV/ result. 
We can motivate the 2/6 spacing of the renormalons by considering the 
connection between the IR renormalons and the absence from the formal 
perturbation theory of non-perturbative vacuum condensates present in the 
OPE. Performing a short distance OPE on the time-ordered product of cur- 
rents in equation (1) one has 

(24) 
where CpT corresponds to the perturbation series of equation (6). There 
are then higher terms consisting of perturbatively calculable coefficient func- 
tions multiplied by non-perturbative vacuum condensates of increasing mass 
dimension. The lowest such is the gluon condensate with dimension 4 and 
a Q~'^ scaling behaviour. There will also be dimension 6, 8, . . . condensates 
corresponding to Q"^, Q~^, ■ ■ ■ behaviour. Crucially, there is no relevant di- 
mension 2 operator and hence no condensate with Q~^ behaviour. These 
condensates are supposedly associated with the IR renormalon singularities 
at z = 4/b,Q/b, . . .. The connection follows from a consideration of the Q^ 
behaviour of the large-order perturbative coefficients dk- Suppose that, for 
fi^ = Q^, dk behaves asymptotically as 

dk^A(-)\\{), (25) 

\Zi / 

where A is a constant and the bracket denotes possible additional powers or 
logarithms of k. Then, assuming a trivial beta- function with one term (c=0), 
the renormalization group fixes the asymptotic behaviour for general /i^. 




■4(^1 (-|«(). (26) 

So such a singularity at z = Zi corresponds to Q~^^^ scaling behaviour. To 
reproduce the Q~^, Q~^, ■ ■ ■ behaviour of the OPE condensates one then re- 
quires Borel plane singularities at z = 4/6, 6/6, . . ., the IR renormalons. 

In QED, one-loop vacuum polarization diagrams with a chain of vac- 
uum polarization bubbles inserted lead to fixed-sign factorial growth and 
UV renormalon singularities ]T^ . These are associated with the Landau pole 



in QED. For the QCD vacuum polarization function one can also consider 
one-loop diagrams with a single gluon line inserted. Applying a cut-off on 
the momentum of this line, inserting the running QCD coupling and inte- 
grating over the high-momentum region, Mueller []TU[ has shown that one can 



explicitly derive the form of the leading UV renormalon at z = —2/6; and, by 
integrating over momenta less than the cut-off, of the leading IR renormalon 
at z = 4/6. 

The final conclusion is that one expects UV renormalon singularities, 
UV£, at z = ze = —2i/b, i = 1,2,3,...; and IR renormalon singularities, 
IR^, at z = Zi = 2i/b, i = 1, 2, 3, . . .. For the specific case of QCD vacuum 
polarization one expects IRi to be absent since, as discussed, there is no 
dimension two condensate. 

These singularities at z = Z£ in the Borel plane should be branch points 



of the form flS 



^ . -4o + ^(l---A,) + 0((l -./-->)') p,j 

where p is a positive integer. The exponent czg, is fixed by the renormaliza- 
tion group in order to give the required Q^ scaling, taking into account a 
realistic two-term beta-function with c 7^ 0. UVi should be the singularity 
closest to the origin (assuming that IRi is absent) and so will give the overall 
asymptotic large-order behaviour of the d^ coefficients. Corresponding to a 
branch point in the Borel plane of the form of equation (27), one has the 
large-order behaviour 

4 = p, 1° , (-) T{k + p + cz,). (28) 

T{p + cze,) \Zi) 

10 



So for UVi we have the leading behaviour 

A. ( h\^ 

dk ^ 







T{p- 2c/ b) 



k''^^k-^^/''k\{l + 0{l/k)) 



(29) 



where we have expanded the second F-function of equation (28) for large k. 
The Ai and higher terms in the numerator Taylor series about z = ze in 
equation (27) are 0{l/k) sub-asymptotic effects. The explicit diagrammatic 
evaluations ||, 0, |T9[ give Aq = |r(2 - 2c/b), p = 2 for UVi, in the MOM 
scheme with ^"^ = Q"^. 

To make contact with the 1/Nf and 1/6 expansions of equations (9) and 
(12) we note that from (8) 



c 
l 



8 62 + 4 5 



11 CyJ^C/p 3 Cyp 

y 62 + 1-j- 



(30) 



Then writing k '^^''^ as e ^cinfc/fe ^^^ expanding we have 



dh 



Ao 



T{p-2c/b) 



-- kP-^k\ 



5 Ca , , 3 Cp , , 

Ink m k 

2 6 2 6 



/25 



7 



+ — ln'A; + -lnA; -# + ••• 






4 



C^ 



62 



We shall define 



An 



T{p- 2c/ b) 



/l02 Ao3 



(31) 



(32) 



where A^r is a sum of multinomials in Ca and Cp of degree r — 1 and Aqi 
is a pure number. One can then obtain the asymptotic behaviour of the 
coefficients to all orders in the 1/Nf and 1/6 expansions. The 1/6 expansion 
corresponds to the successive terms in equation (31). So we have 



d 



di' 



(fc-i) 



Aoi 

Aqi 



_1^ kP'^k\{-^-CA\Tik-^-Cp\nk 



(33) 
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For the l/Nj expansion asymptotics the leading terms will come from the 
binomial expansion of h^ = ( — ^ — -)'', 

where r = 0,1,2,.... Note that an additional factor of k is gained for each 
additional order in the 1/A^/ expansion. The leading term in the large- A^ 
expansion, d[% is given by 

where Aqi is the large- A^ limit of equation (32). Of course Aqi cannot be 
determined from the large- Nf result. 

Assuming a UVi singularity as in equation (27) we therefore apparently 
have the remarkable conclusion that, given an exact large- A^/ (large-6) result 
for df, (dj. ), we can obtain Aqi and p; and hence the asymptotics (up to an 
0{l/k) correction) of the coefficients to all ordersin the 1/A^/ (1/&) expansion 
are determined by equations (33) and (34). We shall check in due course that 
the exact large- A^/ result of equation (18) gives Aqi = 4/9 and p = 2 (MOM 
scheme, /x^ = Q"^) in agreement with other evaluations |^, |^, |TU[. 

The conclusion that the large-Nf result can determine the full asymp- 
totics beyond leading order in 1/A^/ seems too good to be true; and indeed 
recent work by Vainshtein and Zakharov casts doubt on it. These au- 
thors have systematically developed a "UV renormalon calculus" in which 
the leading ultra-violet behaviour of loop diagrams with different numbers of 
chains of vacuum polarization graphs inserted is extracted using an operator 
product expansion. Including one chain they reproduce a UVi result of the 
form of equation (27) in agreement with the result quoted above. Explicitly 
evaluating the two-chain (three-loop) result in a simplified U(l) model they 
find a contribution to the UV renormalon asymptotics which is 1/A^/ down 
on the one-chain result; but has additional powers of k and so dominates 
the one-chain result. Conventional wisdom would have expected additional 
chains to have a suppression by powers of k and hence not to contribute to 
the Aq coefficient in equation (27). Combinatoric factors conspire to modify 
this, however. 
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The UV renormalon calculus results ^ imply that the Borel plane sin- 
gularity at z = —2/6 in equation (27) should be modified to 

^^^' ~ h (1 + 6^/2)7™-2c/. w 

where m labels the number of vacuum polarization chains, m = 1 correspond- 
ing to the previous one-chain results. Counting powers of Nj, one expects 
each additional chain to give a 1/A^/ suppression and hence equation (32) 
should generalise to 

^' -^01+4°^ + ... (37) 



r(7^ - 2c/b) 6—1 

The exponent 7^ will be related to the anomalous dimensions of the operators 
appearing in the operator product expansion of reference 0. There will be 
a number of operators and hence a number of contributions to equation (36) 
for any m. Each such 7^ will have a '^/Nf, 1/6 expansion and we can define 
the leading term pm, typically a positive integer. We select for each m the 7™ 
contribution with largest Pm] and this is the single term displayed in equation 
(36). Providing that Pm+i > Pm + 1 for any m, then (36) and (37) lead to 
the 1/A^/, 1/6 expansion asymptotics 

dt'^ ^ 4r')(-i)'F--^fc! (38) 

with r = 0,1,2,.... The Aqi will consist of multinomials in Ca, Cf of 
degree r. The inequalities on pm are required since, from equation (34), we 
see that, at fixed m, one gains an extra factor of k for each additional l/Nf 
order. The results of (38) imply that one does not get something for nothing 
after all; but requires 0{{l/NfY) results (diagrams with r + 1 chains) to 
obtain the leading asymptotics to this order in the 1/iV/ expansion. 

Having discussed the general expectations for the Borel plane singularity 
structure we return to the exact large- A^j result of equation (18) and exhibit 
what its singularity structure actually is. 



13 



4 Borel Plane Singularities of the Exact Large- 
Nf Result 

Returning to equation (15) we can use Taylor's theorem to write 

*L'^' = ^-(^-2)!(D„-2PH (39) 

where (D„F(x) denotes the coefficient of x" in the expansion of F{x) as a 
power series. 

Pin) = ^^ y ^-'^'' (40) 

"^^"^ 3(2 + .),^,(P-(l + .)T ^'"^ 

which is obtained from equation (16) by writing u = x — 1. Using partial 
fractions on equation (40) and isolating the required coefficient one finds the 
rather simple result 

^L"l 6n-l 6n+l ^, ..^^ f / 1 1 



+ 12E(-l)^+i|7 



(n-2)! 3"+i 6"+i ^^3 

1 2 M / 1 

+ -7^— TTT + 



(£ + 1)2 (£ + 2)2yj V3£ 

-(=|^ + i2g(-i,«{„(^ 

^v(£-i)2 ~ {i-ifnvM) ■ ^^^^ 

Putting in factors of T^N^ = (1/2)"A^? and replacing Nf by (—36), we see 
that successive terms in equation (41) are proportional to (—6/2)", (—6/4)"^, 
(— 6/2£)", (6/4)" and (6/2£)" corresponding exactly to the poles UVi, UV2, 
We, IR2 and IR^ respectively in the Borel plane. There is no term involving 
(l/(— 3))" and so IRi is indeed absent as anticipated. The linear factor 
proportional to n in all but IR2 means that all the poles are double poles 
{p = 2) except for IR2 which is a simple pole {p = 1). Notice that the branch 
point exponent cze is sub-leading in the '^/Nf, 1/6 expansion and so one sees 
poles and not branch points in the large- A^/ limit. 
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Using equation (18) we can then obtain the leading-6 resuh for dk (MS 
scheme, fj,"^ = Q"^): 



f oo 

dfh'^ = k\\j2[Aoii)k + A 

oo 



b ^'^ 

2£. 



(42) 



where 



^o(l) 

Ao{2) 

Ao{e){i>3) 



Boil) 
Boi2) 

Boii) 



^e-/^^i(l) = ye-/^ 




-10/3. 



5i/3. 



2. 
-5£/3. 



(£+1)2 (£ + 2)2^ 

Bi{l)=0; 
0;5i(2)=ei°/3 

-Aoi-i); B,{£) = -A,{-i) {£ > 3) 



(43) 



To obtain the result in the MOM scheme with yU^ = Q^ one simply omits 
the exponential factors in equations (43). The first sum in equation (42) 
generates UV^ singularities of the form of equation (27) and the second sum 
the IR^ singularities. IRi is absent (-Bo(l) = -Bi(l) = 0) as required from the 
OPE; and all poles are double {p = 2) except IR2 for which Bq{2) = giving 
a simple pole [p = 1). For i > 3 there is a curious and unexplained symmetry 
between the residues of UV^ and IR^ with A{£) = —B{—i). As promised the 
coefficient Ao(l) = f^"^^^ is in agreement with other calculations of UVi 
[|, 1^, 0. The coefficient -Bi(2) = e^°/^ is consistent with the result of [|19 
for the first IR renormalon, IR2. 

Since equation (41) is split into contributions from UV and IR renor- 
malons we can sum up these contributions separately. Summing the first 
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three UV terms in (41) one obtains 

^L"i (UV) - 2(C2-2) 2C2-3 



(n-2)!' ' (-3)"-! (-6)" 

4 "-1 



and summing the final two IR terms gives 

^W ,^^, 2(n2-3n + 4-C2) ^(n^ + 3n + 2 - 4(2) 



r(IR) 



(n-2)!' ' (-3)"-! (-6)"-! 

4 '"-1 



■ i ov^ E ^(1 - 2-™)(l - 2"^-")C„.+i. (45) 

The UV and IR pieces separately contain even and odd ^-functions but the 
even (^-functions cancel in the sum of (44) and (45) to reproduce equation 
(17) which contains only odd (^-functions. 

We finally discuss the connection between the Borel plane singularity 
structure of D, which we have discussed extensively, and that of the more 
experimentally-relevant i?, related to it by the dispersion relation of equation 
(4). It is easy to show that in the large-6 limit (c = 0) [|ll| 



BW,) = ^i^^BPK.) (46) 

Since sin(7r6z/2) has single zeros ~ (z — Z() at the same positions as the 
renormalon singularities one finds that renormalon poles of order p in B\D\ 
are converted to poles of order p — 1 in i?[i?]. This implies that the poles 
in -B[-R] are simple poles except for IR2 which was a simple pole in B\D\ 
and hence apparently vanishes [§. The absence of the IRi singularity at 
z = 2/6 in B[D]{z) implies from equation (46) that B[R]{z) must have a 
compensating zero at this position. Brown and Yaffe |jll[ considered this 
unlikely and hence cast doubt on the absence of IRi. The exact large- iVj 
result shows that there is indeed no IRi singularity in B[D] and hence such 
a zero is present in B [R] . 

The leading asymptotics of the coefficients r^ will be given by UVi for 
B[R]{z). Expanding around z = —2/6 we have 



^^^K-l)-«((-l)l '-) 



7ibz/2 
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and hence from equation (46) we find that the asymptotic behaviour of r^ 
is given by changing p — > p — 1 in the results for dk (equations (33) and 
(34)). Even assuming the more comphcated UVi structure of equation (36) 
one simply changes Pm —^ Pm " 1- This implies that on very general grounds 
one expects 

g^^(l + 0(lA)) (48) 

so that the r^ coefficients grow more slowly asymptotically. 

We conclude by noting that exact large- A^j results do exist also for other 
QCD observables. In particular, reference |^ contains a leading-A^j result 
for the perturbative coefficients of the radiative corrections to the Gross- 
Llewellyn-Smith (GLS) sum rule. This reveals that in the Borel plane the 
GLS corrections have simple poles at z = ±2/6, ±4/6. So the first two UV 
and IR renormalons are present but the remaining renormalons are absent 
to leading order in 1/Nf. It would be interesting to compare this with the 
OPE structure for these corrections. 



5 Vanishing IR Renormalons 

The singularity structure in the Borel plane for SU(A^) QCD will change as N 
and Nf are varied. In particular, ii Nf approaches llA^/2 from below then b = 
{UN — 2Nf)/6 will approach zero from above; and, as more flavours of quark 
are added, the IR^ and UV^ singularities, which are spaced at intervals of 2/6, 
will move outwards away from the origin in the 2;-plane. The instanton/anti- 
instanton (iT) singularities remain fixed at z = 4, 8, 12, . . . independent of A^ 
and A'^^. When 6 = 1/2 {Nf = 15 for SU(3)), IRi will be at the same position 
as the leading IT singularity at z = 4; and for fiavour saturation (maximum 
Nf for which 6 > 0, Nf = 16 for SU(3)) 6 = 1/6 and the leading singularity 
on the positive axis will be the iT at 2 = 4. 

For SU(3) QCD a remarkable phenomenon first occurs at 6 = 1/2 (A'^^ = 
15). One finds that the branch point exponent of IR^, 2ci/b in equation (27), 
becomes a negative integer and the structure of IR^ in the Borel plane is then 

/ ^ N 88i-p 

B[D]{z) = Ao[l - -) . (49) 

The IR^ singularity disappears provided that p < 88i. For the particular 
case of the D-function we have apparently p = 1 or 2 and so all of the IR 
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renormalons disappear at 6 = 1/2 and the only singularities on the positive 
z-axis are those due to instantons. Notice that the UV renormalons are still 
present but become poles. 



As first noted by White |12, 13| the IR renormalon singularities also 



disappear for Nf = 16 flavour-saturated SU(3). The exponent 2ci/b again 
becomes a negative integer and the IR^ structure in the Borel plane is 

B[D]{z) = Ao[l - —) . (50) 

So for p < 906i the IR^ singularities disappear. 

The implication is that Nf = 15 and 16 SU(3) QCD are very special 
instanton-dominated theories. At precisely the point when instantons be- 
come the leading singularities all the other IR renormalon singularities on 
the positive z-axis vanish. It is interesting to ask if this scenario is unique 
to SU(3) or can be realised for other values of A^. 6=1/2 will occur for 
integer Nf only for A^ odd. For A^ even the iT singularity becoming leading 
and flavour saturation are telescoped into the single value b = 1/3. 

Considering A^ odd first, we require that c/b is integer for b = 1/2; 



25N^ + 13N^ + UN - 3) 

■ (51) 



b=i 4:N 

" 2 



A necessary condition for this to be an integer is that 3 = (mod A^) , which 
uniquely fixes A^ = 3. Then (c/6)|^^i = —44. 
For A^ odd and flavour saturation b = 1/6; 



225A^=' + 39A^2 ^ ggjy _ 9) 

. (52) 



=1 4A^ 

6 



A necessary condition for this to be an integer is that 9 = (mod A^) so 
A^ = 3 or 9. The A^ = 9 case gives (c/6)|b=i = -4444; and, for A^ = 3, 
(c/6)|,= i=-453. 

6 

For A^ even and 6=1/3 



(-225A^3 ^ 78JY2 ^ ggjy _ ^g) 

.a = m • (^^) 



"3 



A necessary condition for this to be an integer is that 18 = (mod N) and 
so iV = 2, 6 or 18. A^ = 6 is the only case for which c/b is an integer and 
then {c/b)l^i = -1884. 

So for flavour-saturated SU(A^) IR renormalons are absent only for A^ = 
3, 6 or 9. For SU(9) the b = 1/2 case where the I T singularity becomes leading 
still has IR renormalons. 

We conclude that SU(3) QCD is a very special theory from yet another 
point of view. 

6 Conclusions 

In this paper we have discussed and sought to extend our knowledge of 
the Borel plane singularity structure of the Adler D-function (QCD vacuum 
polarization). This singularity structure succinctly encodes the large-order 
asymptotic behaviour of the perturbation theory coefficients. 

We pointed out that an expansion of the perturbative coefficients in pow- 
ers of b, the ffist QCD beta-function coefficient, rather than in Nf, was 
natural when comparing with QCD renormalon expectations. We further 
noted that, if the leading UV renormalon indeed has the expected structure 
of a simple branch point, then knowledge of the perturbative coefficients to 
leading order in 1/Nf, 1/6 allows the large-order behaviour to all-orders in 
'^/Nf, 1/6 to be inferred, equations (33) and (34). This seems unlikely and a 
more complicated structure was proposed, equation (36), which is consistent 
with the UV renormalon calculus results of Vainshtein and Zakharov [^. 

Using the exact large-A^/ result for the D-function of reference [^j we 
exhibited the explicit singularity structure in the z-plane and found the ex- 
pected UV and IR renormalon singularities. They appear as poles in the 
large- A^/ limit. In particular, the ffist IR renormalon, which would corre- 
spond to Q~^ behaviour not present in the OPE, is absent. We gave explicit 
expressions for the residues at all of these poles (equations (43)) and unex- 
pectedly found those for the third and higher UV and IR renormalons to be 
symmetrically related. We were also able to sum up separately the UV and 
IR renormalon contributions in closed form (equations (44) and (45)) and 
obtained expressions containing even and odd (^-functions. When UV and 
IR contributions are combined the even (^-functions cancel. 

We finally noted that in flavour-saturated SU(A^) QCD the IR renor- 
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malons are absent for A^ = 3, 6 and 9. These theories then have ambiguities 
dominated by instantons. For SU(3) the IR renormalons first disappear when 
Nf = 15 (6 = 1/2), at which point the IT singularity is in the same position 
as IRi and becomes leading. 
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